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Abstract 

Self-phase modulation of spherical gravitational wavepackets propagating in a flat space-time in 
the presence of a tenuous distribution of matter is considered. Analogies with respect to similar 
effects in nonlinear optics are explored. Self phase modulation of waves emitted from a single 
source can eventually lead to an efficient energy dilution of the gravitational wave energy over an 
increasingly large spectral range. An explicit criterium for the occurrence of a significant spectral 
energy dilution is established. 
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I. INTRODUCTION 



It is well known that the Einstein's equation describing gravitational waves is strongly 
nonlinear Nonlinear wave processes similar to those observed in Optics can then 

eventually occur. Recently, the possible occurrence of self-phase modulation, harmonic 
generation and nonlinear wave mixing was considered 13,01 . 

This is an important issue in two different aspects. First, from a theoretical point of 
view it is important to explore and to understand the similarities and differences between 
gravitational wave phenomena and effects in Nonlinear Optics that have been tested in 
laboratory. Second, and in more practical terms, this is important in what concerns the 
possible detection of gravitational waves. Detectors have been designed and built under 
the assumption that the frequency spectrum of gravitational waves emitted by astronomical 
objects is conserved, and that the wave intensities decrease as the inverse of the square 
of the distance. However, the nonlinear processes can eventually lead to spectral energy 
dilution such that the energy density received within the detectable frequency bandwidth is 
significantly decreased. Other processes eventually contributing to energy dilution could be 
due to the coupling with plasma waves 0] and with photons p, Q| • 

For waves emitted from a single astronomical source, the main nonlinear effects that 
can occur are self-phase modulation (for short pulses, of the order of a few cycles) and 
harmonic cascades (for longer pulses). Attention was however called to the fact that, for 
parallel propagation, the strong nonlinearities associated with empty and flat space-time 
exactly cancel each other Q]. The existence of such a negative result is apparently due 
to the absence of gravitational wave dispersion. Note, however, that anti-parallel wave 

are not relevant to waves 

emitted by single sources. 

In this work, we return to the problem of parallel wave interactions. We will focus on 
self-phase modulation of spherical waves emitted by isolated sources. For simplicity, we will 
consider a fiat space-time, filled with a tenuous distribution of matter j^. The matter dis- 
tribution will guarantee the existence of wave dispersion, which is an important ingredient 
of self-phase modulation in Nonlinear Optics The nonlinear wave equation is estab- 

lished in Section II, and the dispersion properties of linear waves is discussed in Section III. 
Nonlinear evolution of a spherical gravitational wavepacket is studied in Section IV, where 
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a necessary criterium for the occurrence of a significant amount of self-phase modulation is 
estabhshed. Finally, in Section V we state the conclusions. 



II. NONLINEAR WAVE EQUATION 

We consider propagation of small amplitude gravitational waves, in a region of space-time 
where we have a tenuous distribution of matter j^. We can then, in a first approximation, 
neglect the background field curvature. 

We are considering a flat space-time, perturbed by a small amplitude gravitational wave. 
This can be described by the metric tensor elements 

9ij = Vij + hij, (1) 

where \hij\ <^ 1 represent the gravitational wave, and rjij are the metric tensor elements of 
fiat space-time 

?7oo = 1 , Vii = -1 = 1,2,3) , r]ij = {i ^ j). (2) 
In this case, we can derive from Einstein's equation the following nonlinear wave equation 

n^h,k = -2KS,,+2Rfl\ (3) 

where k = {8tt/c'^)G, and G is the gravitational constant, and is the d'Alembert operator 
in the usual form 

□2 = djg. = r]'^didj. (4) 

where we use di = d/dx^. The nonlinear term Rf^ contains third order nonlinearities of the 
Ricci tensor. It can be explicitly written as 

— {dkhpi + dihpk — dphik) dihnp] ■ (5) 

(2) 

Here we have neglected the second-order nonlinear term, R^j^ = 0. Note that the second 
order part of the Ricci tensor does not vanish identically unless the response due to the 
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pseudo energy- momentum tensor is taken into account 4". Combined with the original 
perturbation, these second order terms contribute to the self phase modulation of the wave. 
The equation studied here may therefore be regarded as a model equation for the nonlinear 
dynamics of gravitational waves. In equation we have also included the linear dispersion 
term associated with the matter distribution uSij, where Sij is related with the energy- 
momentum tensor Tij by 

Sij = Tij - ^riijT (6) 

where T = T? is the trace. 

III. LINEAR DISPERSION RELATION 

Let us first consider the properties of a linear wave, propagating radially from a given 
point source. In order to discuss the dispersion properties of this wave we first consider the 
linearized wave equation 

O^hik = -uSij, (7) 
and assume a plane wave solution of the form 

hij = eijAexplignx"^], (8) 

where i = v^—T, A is the amplitude, eij is unit polarization tensor such that e*jeij = 1, and 
Qn are the components of the four-wavevector. If the scale of variation of the amplitude A 
is much larger than the typical wavelength, we can use djhik = iqjhik, and write the linear 
wave equation as 

il^'^qjQnhik = 2KSik. (9) 

The perturbed energy-momentum tensor can be considered proportional to the local 
amplitude of the gravitational wave: Sik = WikA, where the tensor Wik depends on the 
properties of the medium. We are then lead to the following linear dispersion relation 



v^^'qjqn = w, (10) 
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where we have used w = 2ne^^*Wik- Particular examples of w can be found in the literature. 
For instance, the cases of a cold dust cloud |3l and of a magnetized plasma are well 
established and don't need to be explicitly given here. The contribution to dispersion from 
the background curvature has also been calculated for various space-times, see e.g. 

In order to deal with spherical waves from a localized source it is appropriate to use a 
spherical coordinate system {r,6,(f)), such that: 

x° = ct, = r, = e, = 0. (11) 

For waves propagating in the radial direction we have to replace solution (jSJ by the 
following wave solution 

hik = (^ik^ exp{iqoX° + iqix^) = tik^ exp{iqr - iQt), (12) 

where a is the new amplitude, q = qi and = —qoc. Because the nondiagonal components 
7]^'', with i ^ k, are equal to zero, we can easily transform the dispersion relation (fTUI) into 
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-q'] =w{r,e,<f)), (13) 



where we have retained the possibility of a non-uniform distribution of matter. We can 
see that the matter distribution can change the phase velocity of the gravitational wave, 
according to 



For the group velocity, we have 



dVt qc 

It is clear that, even if w is independent of q, the presence of a small amount of matter 
leads to a group velocity wave dispersion. It is known from Nonlinear Optics that wave 
dispersion is an essential ingredient of self-phase modulation It is the absence of 

dispersion that eventually explains the otherwise counter-intuitive result that self-phase 
modulation is absent for plane gravitational waves propagating in empty fiat space-time 
The inclusion of matter is thus an essential ingredient of the present study. 
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IV. NONLINEAR WAVE PROPAGATION 



We can now examine the possibility of a given radial wave, satisfying the above linear 
dispersion relation, to interact with itself, due to the nonlinear contributions contained in the 

(3) 

term R^j^ ■ The nonlinear contributions of Sik could equally be included, but for simplicity 
they are neglected here. The existence of nonlinear wave coupling implies that the wave 
amplitude a in equation ()12|) can no longer be a constant, and is replaced by a slowly 
varying function of r and t. This means that we now have 

djhik = {iqj + -dja)hik. (16) 
a 

So, we can write 



n^hik ~ [-V^^'qjqn + In a]hik. (17) 

Assuming that the above linear dispersion relation still holds, we can cancel the first of 
these terms with the linear contribution from Sik- For the second term, we can write it as 

W'^Qjdn = -iqdr - i^dt = -iq{dr + -dt) (18) 

where v = c^/vj is the group velocity. 

In order to establish the nonlinear equation for the slowly varying amplitude a we now 
use an approximate expression for R-f^ , where only the terms oscillating at the frequency of 
the wave Q are retained 



R^-^i: -^-<l')^-eM^1r-^m■ (19) 



1 /n'^ _ ^\ lapa 

2 ^ / r 

In deriving this expression we have added to the solution (fT^ its complex conjugate, in 
order to adequately describe a real wavepacket. For a given line of sigth between the source 
at r ~ and the eventual observer at a finite distance r, we can use equation (|T^ with fixed 
values of 9 and 0. Replacing in the nonlinear wave equation (jH} we obtain 



1 „ X , Aa 
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iq(dr -\ — dt)a = w(r)^r-a (20) 

V r'^ 

This nonlinear equation for the slow wave amplitude clearly shows that the nonlinear 
effects disappear in the absence of matter, w{r) = 0, as noticed previously j^. Let us make 



a variable transformation from the pair (r, t) to {z, r), where we define z = r — vt and t = t. 
We have then dr = dz and dt = dr — vdz- Replacing this in the above equation, we get 



.w{z,t) lap 

d^a = -I ^ ' ^ V , a 21) 

q r^[z,Tj 

This equation is satisfied by a solution of the from 

a{z, r) = a{z) exp[i(f){z, r)] (22) 
with the phase function determined by 



""^ w(z, t') la 

[Z.T] =0a- I 



2 



-dr' (23) 



q r'^[z,r 

This solution represents a gravitational wavepacket propagating spherically with a nearly 
constant envelope a{z) and a variable nonlinear phase. The wave frequency shift AQ will 
be given by the derivative of this phase with respect to the time variable t 'll|. Neglecting 
the small variation of the distance r{z,T) and the matter distribution w{z,t) inside the 
wavepacket envelope, this means that AQ will be essentially due to the variation of the 
energy distribution |a(2;)p with respect to time. But this envelope is only a function of 
z = r — vt and we can use 



dM^)\^ = -^dM^)\^ (24) 



We can then state that 



An{T) = -vdz(l){z,T). (25) 

Noting that, for short wavepackets, the variation of the matter dispersion term w{z, r) 
and distance with respect to the source r{z, r) can be negligible, we replace them in the 
expression of the phase by their central values w{t) = w{z = 0, r) and r(r) = r{z = 0,r). 
This leads to the following expression of the frequency shift occurring inside the wavepacket 
envelope 

mr) = f ^^^dM^)\'dr' (26) 
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Here we notice that the distance travelled by the wavepacket can be written as r(r) = 
v{T')dT' CT. Neglecting the possible slow change on the shape of the envelope over 
distance, we can finally write the above expression as 

Anir) -dMz)\' r^dr' (27) 

q Jo ^ 

In order to understand the physical meaning of this result, let us consider the simple case 
of a uniform distribution of matter along the entire line of sight. We can use w{t) ~ = 
const., and get for the frequency shift, after a distance r c:^ ct travelled by the wakepacket 

AQ{r) ~ --wod,\A{z)W (28) 
1 

where we have used the local spherical wave amplitude A{z) = a{z)/r, observed at a distance 
r from the source. This linear dependence of the frequency shift with time, or with the 
travelled distance, was found previously for linear propagation and is well known from 
Nonlinear Optics Q|. Here, however, the frequency shift is proportional to the square of 
the local amplitude, which means that this effect can only the significant if it occurs over 
short distances, not far away from the emitter. This feature is specific of spherical waves 
propagating in uniform media. 

Another interesting case is that of a non-uniform matter distribution where the 
wavepacket propagates accross a succession of localized clouds, at distances (with 
i = 1...N) from the source, and widths Ar^ ^ r^. We can then transform equation pTj) in 

AQ{r) ~ -- Vw;A|A,(^)|2Ar, (29) 
q ^-^ 

where the local envelope amplitudes are determined by Aii^z^ = a{z)/ri. Again, the strongest 
contribution to the total frequency shift will result from the clouds located nearest to the 
source, supposing that they all have similar matter densities. 

Let us assume that the source emits a gravitational wavepacket with n cycles. Its width 
with be 6z ^ 27m/q. And the maximum frequency shift associated with the closest cloud 
(i = 1) will be of the order of 

AQmax ^ -wi-^ 30 

q dz 
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We can also write wi = {Q/c)'^a, where a <^ 1 is a small dimensionless factor. Noting 
that Q ~ qc, this allows us to establish the necessary condition for a large frequency shift, 
leading to a significant spectral energy dilution, as Afl > fl, or equivalently 

a\A,\'^ = ^\A\\qAn)>l (31) 

Notice here that (gAri/27r) is the number of wavelengths over the cloud width. This can 
be a very large number. For instance, for Q = 10'^ (wave in the KHz range) and a width 
of Ari = lO^^parsec, we get (gAri) = 10"'^'^. This means that, for a short pulse {n < 10), 
dense enough cloud (a > 10"'^) and close enough source {\Ai\ > 10~^), the above criterium 
could be satisfied. This confirms, in more solid grounds, the suggestion previously made 
that self-phase modulation could eventually take place. 

Some further discussion on the size of the effect of self-phase modulation may be in place. 
In principle, close to the source the self-phase modulation could be stronger than the value 
presented above, not only because of larger gravitational wave amplitude, but mainly because 
of the significant angular dependence of the metric. This can be seen as follows. Observe 
that terms proportional to — c^q^ in Eq. (|19|) implicitly contains an angular dependence, 
due to the fact that for spherical waves q = q{r) is found from the spherical Bessel functions. 
Separating variables in spherical coordinates, it can be inferred that Q'^—c'^q'^ ~ £(£— l)c^/r^, 
i being the mode number of the associated Legendre polynomial. Thus, for high mode 
numbers, the nonlinear modification may increase by a substantial amount. 



V. CONCLUSIONS 



Nonlinear wave propagation of spherical gravitational waves was considered in this work. 
The possible occurrence of self-phase modulation was discussed. The case of short wavepack- 
ets emitted from a point source in fiat space-time was examined, where a tenuous distribution 
of matter was retained in order to guarantee linear wave dispersion, which is a necessary 
condition for self-phase modulation to occur. An explicit criterium for a significant spectral 
energy dilution due to self-phase modulation was established. It leads to the conclusion 
that the occurrence of self-phase modulation due to matter distribution very close to the 
gravitational wave source is plausible. 

In contrast with what could occur with plane wave propagation, for spherical waves the 
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contributions of phase modulation over distance decay very rapidly with distance from the 
source, due to wave amplitude decrease. For this reason, self-phase modulation is dominantly 
occurring very close to the emitter. Notice however that the region where this effect takes 
place is not necessarily the region where it can be observed, because the expanded wave 
spectrum will then propagate far away without further changes. 

The efficiency of the self-phase modulation process is directly dependent on wave dis- 
persion, which is a consequence of matter distribution. Curvature of space-time would also 
contribute to wave dispersion and would enhance the process. If, instead of spherical wave 
emission we have some kind of directionality, the wave amplitude decay will be smaller and 
phase modulation will also increase. Another source of nonlinearity is the energy-momentum 
tensor, or the matter distribution itself, which was not retained here. Space-time curvature, 
directionality effects and energy-momentum nonlinearities will eventually lead to more fa- 
vorable criteria for the occurrence of self-phase modulation of gravitational waves, and will 
be considered in a future work. 
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